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A simple example is used to show that renormalization group limit cycles of effective quantum 
theories can be studied in a new way. The method is based on the similarity renormalization 
group procedure for Hamiltonians. The example contains a logarithmic ultraviolet divergence that 
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the limit cycle with an entire hierarchy of bound states. 
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I. INTRODUCTION 

A model Q invented in the front form of Hamilto- 
nian dynamics @j|3| has been used in its generic matrix 
version to argue [4| that cyclic dependence of coupling 
constants on cutoffs in renormalization group (RG) pro- 
cedure :5| may actually be commonplace for quantum 
mechanical systems whose Hamiltonians require renor- 
malization. In independent studies, many examples of 
physical few-body systems with short range interactions 
and large scattering length have been considered in which 
a cycle structure may occur with universal features 0. 
In the theory of three-body systems in nuclear Q and 
atomic 8] physics, the limit cycle is associated with the 
existence of a sequence of bound states with binding en- 
ergies forming a geometric series. The cycle has also 
been discussed in the context of the Bose-Einstein con- 
densation yj and superconductivity [lfj. Even the ele- 
mentary non-relativistic Hamiltonian for one particle in 
the potential ~ 1/r 2 , when properly regulated, exhibits 
the cyclic behavior [Tl|. Concerning particle theory, the 
cyclic behavior of effective nuclear interactions could oc- 
cur if the masses of 7r-mesons, which influence the shape 
of nuclear potentials, were increased by about one third. 
Since the masses of 7r-mesons can be linearly related to 
the masses of up and down quarks using the Gell-Mann- 
Oakes-Renner relation, the nuclear theory suggests that 
the masses of light quarks in the Standard Model are 
near the range where QCD may develop an infrared (IR) 
limit cycle 0> El • I n order to determine if the IR limit 
cycle does indeed occur in QCD, the binding mechanism 
for quarks and gluons will have to be understood much 
more precisely than it is so far. While the mechanism 
of binding in just three-body systems has fascinated re- 
searchers for many years |l4l and a large body of lit- 
erature exists on the subject studies of the RG limit 
cycle appear to be in their infancy [T^ . Quantum three- 
body calculations that trace the dependence of coupling- 
constants on cutoffs tend to be complex ^^|. Thus, the 
simple matrix example provides useful insights con- 
cerning universality in Hamiltonian quantum mechanics 
with RG limit cycle pf . 

So far, all studies of limit cycles known to the author 



have been carried out in one of two ways. The first way, 
reminiscent of the concept of charge renormalization in 
quantum electrodynamics [20|, is to regulate a diverg- 
ing model with a cutoff and solve for some observable, 
such as a scattering amplitude or a bound-state energy, 
which is postulated not to depend on the cutoff. When 
the cutoff parameter is varied, it is possible to keep the 
observable fixed if instead a coupling constant is allowed 
to vary (there may be more than one constant to vary). 
Limit cycle means that the coupling constant is a periodic 
function of the cutoff parameter when the latter increases 
to infinity. Most of the quoted literature discusses exam- 
ples of a limit cycle in this way. But if it is known that a 
theory being considered is only applicable to phenomena 
of a limited range of scales, the limit of infinite cutoff 
becomes logically questionable. The limit may not exist 
if the required coupling constants diverge in or on the 
way to the limit. The question of how to guarantee that 
all observables of interest are independent of a changing 
cutoff when it is finite (instead of being infinite) is hard 
to answer proceeding along the first way (unless one deals 
with a model in which the introduction of a small number 
of finite cutoff-dependent coupling constants is sufficient 
to completely eliminate the cutoff from equations that 
describe the model). 

The second way is to apply Wilson's renormalization 
group procedure to a diverging model with a cutoff and 
compute the number of required constants and their 
shapes as functions of the cutoff 0, . The required 
number of constants may be very large, or infinite, which 
is a typical situation for finite cutoffs (see below). The 
calculation involves a reduction of the number of degrees 
of freedom that are explicitly kept in the dynamics, and 
an evaluation of the change of the dynamics due to the 
reduction. Precise calculations can be done using dis- 
crete degrees of freedom and the number of the relevant 
degrees of freedom can be kept finite. The boundary of 
the set of the retained degrees of freedom plays the role 
of a new finite cutoff that is much smaller than the cut- 
off of the initial rcgularization (in the case of fixed-point 
criticality, finite cutoffs only mean that the system can 
be described by the same model at all scales). The fi- 
nite cutoff defines the domain of an effective theory. If 
the calculation is precise, then by construction the finite 
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boundary cannot have any effect on the solution and can 
be set arbitrarily. However, the finite boundary limits the 
effective theory to a small space where the only source of 
dependence on the initial rcgularization cutoff is the ex- 
plicit, computable dependence of the matrix elements of 
the effective-theory Hamiltonian on the initial cutoff. Us- 
ing this information, it is possible to identify the changes 
in the setup of the initial regulated model that guarantee 
that the Hamiltonian matrix elements in effective theo- 
ries with finite cutoffs do not depend on the initial cutoff. 
All predictions of such effective theory are then indepen- 
dent of the initial cutoff. It may happen that the presence 
of new terms, called counterterms, amounts to a replace- 
ment of the initial coupling constants with the new ones 
that depend on the initial cutoff. In principle, the out- 
come of the procedure is a family of effective theories 
with finite cutoffs and predictions that are guaranteed to 
be independent of the finite cutoffs. If a complete the- 
ory must contain many constants depending on the finite 
cutoff, additional conditions, such as the requirement of 
symmetries |23j . may lead to finite correlations among 
otherwise unconstrained finite constants in the effective 
theories. One can consider the limit of an infinite num- 
ber of degrees of freedom by simultaneously increasing 
the initial cutoff and using the RG transformation over 
an increasingly large range of scales to obtain effective 
theories with the same size of finite cutoffs. The limit 
cycle is then the ultimate RG orbit on which the effec- 
tive coupling constants move indefinitely in the case of 
the infinite number of initial degrees of freedom. To ap- 
proximate the limit cycle with increasing precision when 
the number of full cycles in the RG procedure increases, 
the initial regulated model must be set up in such a way 
that all unwanted terms that cause departures from or 
impede approach to the limit cycle in the RG procedure 
are eliminated [19(. 

This article concerns a new way to describe the cyclic 
behavior of effective theories, closely related to but dif- 
ferent from the two ways mentioned above. Namely, the 
Hamiltonian matrices of effective theories are calculated 
using a renormalization group procedure called similar- 
ity |24l |25| . In the new procedure, no states are elim- 
inated from the effective dynamics. Instead, the basis 
states are changed in such a way that the effective inter- 
action Hamiltonian can dynamically couple only those 
basis states whose kinetic energies differ by less than an 
arbitrarily prescribed scale A. This scale plays the role 
of a new renormalization group parameter. A varies from 
infinity in the initial regulated model down to a finite 
value in an effective theory. When the initial regular- 
ization is being relaxed, the Hamiltonian matrix eval- 
uated in the effective basis corresponding to a finite A 
can lead to a dependence of perturbatively calculated 
observables on the initial regularization only through the 
explicit dependence of the Hamiltonian matrix elements 
on the initial cutoff. This explicit dependence is used to 
find the counterterms in the Hamiltonian of the initial 
model. Once the counterterms are found, the effective 



theories are parametrized by finite values of A and all 
have the same predictions that are independent of the 
initial regularization. The dependence of the effective 
Hamiltonian matrices obtained using the similarity RG 
procedure on A is described in the following sections in 
the case of the same generic example that was previously 
studied in Refs. 3 and using other methods. 

The similarity procedure requires a generator (see next 
section) and details of the cyclic behavior of effective 
theories depend on how the generator is chosen. Since 
the main features of the cycle that are related to the 
spectrum of a renormalized model are not sensitive to 
the choice of the generator, variations in the cyclic be- 
havior due to variations in the similarity generator will 
not be discussed in any detail. The explicit discussion 
in this article is focused on the case of a very elegant 
generator taken from Wegner's flow equation for (par- 
tial, in the case of a degenerat e sp ectrum) diagonaliza- 
tion of Hamiltonian matrices [261 l27j . Similarity RG 
procedures using other generators may have better con- 
vergence properties in perturbative calculations |28| , but 
this is of no significance to a qualitative understanding of 
the limit cycle whose description here is based on a non- 
perturbative numerical solution. The model considered 
here does not have degenerate eigenvalues and Wegner's 
generator leads to complete diagonalization when the RG 
procedure is carried out down to energy scales below the 
smallest eigenvalue. 

Wegner [2^| has recently pointed out that differential 
equations similar to his had been considered earlier for 
various purposes (other than a RG procedure) by mathe- 
maticians |3(J,|3l],|32]. The stud y de scribed here is limited 
to the context of RG procedure [2J,|25j . The study shows 
that the limit cycle behavior of effective theories resem- 
bles in part the behavior of asymptotically free ones. 
Since the similarity approach can be applied to gauge the- 
ories (e.g., see [33 and later literature on the light-front 
formulation of QCD), the resemblance observed here in 
the simple model is also of interest in particle physics. 

The model studied here is so simple that the cyclic de- 
pendence on A of a large number of matrix elements of 
effective Hamiltonians can be approximately described in 
terms of just one function of A. This function is called 
the renormalized coupling constant, g\, because it corre- 
sponds to and is closely related to renormalized coupling 
constants in effective Hamiltonians of quantum field the- 
ories. The main result is that each and every RG cycle 
consists of a sequence of gradual changes of g\, like in 
asymptotically free theories, followed by a sudden large 
change that brings the coupling constant back to its ini- 
tial value in the cycle. In principle, this pattern repeats 
itself indefinitely when A changes indefinitely. Each new 
cycle corresponds to a new bound state whose binding 
energy is on the order of A at which the sudden change 
occurs. 

The paper is organized as follows. Section ITT1 describes 
the model and produces the cyclic behavior of g\ . Section 
IIIII explains some elements of the mechanism of change 
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that g\ exhibits in every cycle. Section IIVI concludes 
the paper with comments concerning implications of the 
model study. 

II. MODEL 

The model Hamiltonian is constructed for the purpose 
of qualitative understanding of the logarithmic similar- 
ity RG effects in the presence of bound states, which is a 
typical situation in the case of RG limit cycles. In order 
to construct a model that can be solved precisely, each 
and every energy scale is represented by just one state. 
The energy scale is defined by the spectrum of a Hamil- 
tonian Hq whose eigenvalues are assumed to be given by 
the formula E n — b n , where b > 1 and n is an integer. 
This means that the spectrum of Hq is equally spaced on 
the logarithmic scale with step In b (the unit of energy 
is set to 1). The corresponding eigenstates, denoted by 
|n), form an orthonormal basis, {m\n} = 5 mn . The in- 
teraction Hamiltonian in the model, Hj , is defined by its 
matrix elements in this basis. 

Building models with such discrete basis is not a new 
idea 0, [231 ■ ^ is also a starting point for consideration 
of the case of b ^> 1. The case of large b is not accurate 
physically but it can further simplify the mathematics. 
There may exist complex physical systems where one can 
gain new insights by considering the case of very large b. 
For example, such strategy was instructive in the case of 
the Kondo problem |34j . although quantitative numeri- 
cal results were obtained for b around 2. There is hope 
that some way can be found to make use of a large b in 
studies of light front theory [3j|. There, one also needs 
to find ways to drastically reduce the number of discrete 
degrees of freedom in dimensions other than the energy 
dimension. 

The initial model Hamiltonian, H = Hq + Hj, is de- 
fined by its matrix elements ^| 

H mn (g, h) = (E m E n ) 1/2 [S mn - g - ihs mn ] . (1) 

For m = n, 6 mn = 1 and s mn — 0. For m ^ n, 5 mn = 
and s mn = (m— n)/\m— n\. The interaction Hamiltonian, 
Hi = H — Hq, vanishes when g = h = 0. 

The Hamiltonian of Eq. has a generic ultraviolet 
structure with a logarithmic divergence to which both 
real and imaginary parts of the matrix elements con- 
tribute. Initial regularization is imposed by the cutoff 
E n < A. (In the same cutoff was denoted by A. The 
notation is changed to avoid confusion with A to be in- 
troduced below.) If one sets A = b N , the cutoff means 
n < N and a logarithmic divergence in A means a linear 
divergence in N. When h = 0, the perturbative algebraic 
similarity procedure [24[ produces interaction Hamilto- 
nian with matrix elements of the form g\ (EmEn) 1 / 2 . 
This means that for h = the counterterm that removes 
the dependence on A — > oo from the model eigenvalues 
is obtained by replacing the constant g with a cutoff- 
dependent coefficient g&. In fact, g& exhibits asymptotic 



freedom as a function of A when h = 0. When h ^ 0, gA 
exhibits a limit cycle behavior (or chaos) 0. h does not 
depend on A. 

The similarity procedure yields a family of effective 
theories with Hamiltonian matrices labeled by A. The 
family is described by a solution to the equation of the 
form 

^H X = [F X {H X },H X ] , (2) 

where F\{H\} denotes the generator of the similarity 
transformation. The initial condition is set at A = oo (in 
a numerical study, one can work with inverse of A and 
start from 0, or the infinity is replaced by any number 
much greater than A). In the model studied here, H^ is 
given by H in Eq. with g replaced by gA- The case 
with h — has been originally studied using Wegner's 
form of the generator F\{H\} in [3^, |3||. Writing the 
generator F\{H\} in terms of its matrix elements as 

(m\F x {H x }\n) = f mn (D m - D n )H Xmn , (3) 

where D m — H\ mm is the diagonal matrix element num- 
ber m of H\ , one obtains the following set of equations 
for all matrix elements of H\, 

d 

d\ Hxmn - (4) 

N 

^ [fmk(D m — Dk) + fnk{D n — Dk)]H\ m kH\kn ■ 

k=M 

The lower bound on the summation range, M , is an ad- 
ditional parameter introduced here to simplify numer- 
ical analysis of the low energy region. It sets an in- 
frared bound on the lowest energy included in the model, 
E M = b M . The lower bound is very small for large neg- 
ative values of M. The actual size of M is irrelevant to 
almost all aspects of the limit cycle behavior discussed 
here. The only exception are practical numerical issues 
and the range of A over which one can easily observe the 
cycle using personal computers. M was set to —25 in all 
numerical examples described below. 

Wegner's flow equation j2(| is obtained for f mn = 
dl/dX = — 2/A 3 . This choice amounts to the change of 
the RG variable A to Wegner's parameter I = 1/A 2 . 

The behavior of (N — M + l) 2 matrix elements of H\ 
as functions of A results from solving the coupled nonlin- 
ear differential Eqs. (@J. Numerical experiments showed 
that in order to clearly discern logarithmic effects for 
small coupling constants <?a and h (on the order of 0.1 or 
smaller) using a contemporary personal computer for less 
than an hour, one needs about 30 basis states with b on 
the order of 2. This gives about 1000 functions to trace. 
This number is sufficient for finding dominant features 
in the similarity RG evolution of the coupling constant 
because the behavior of the entire solution can be quali- 
tatively understood by writing (f mn disappears from the 
formula in the case of Wegner's equation) 

H\ mn = \J E m E n A mn exp [— (E m — _E„) 2 /A 2 ] , (5) 
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with matrix A mn of the form 

Amn $mn 9\ ^mn ■ (p) 

The only term depending on A is the coupling constant 
gy, h remains constant. 

Equation © is not precise. For example, using pertur- 
bation theory in the coupling constant <?a, one can find 
other terms in A mn , such as — 2g\ (E m —E n ) 2 /\ 2 . In fact , 
one also finds other functions of the subscripts m and n 
and the width parameter A, and none of these structures 
is present in Eq. Numerical calculation shows (see 
below) that these structures combine to produce the be- 
havior of matrix elements of H\ that is qualitatively cap- 
tured by Eqs. (JSJ) and ijfJJl. All matrix elements can be 
described using variables q mn = [E m — E n )/\. When A 
decreases, q mn for fixed m and n increases and the corre- 
sponding matrix element falls off. The dominant fall-off 
pattern is approximately reproduced by the simple Gaus- 
sian function that multiplies A mn in Eq. (jSJ. Although 
for matrix elements that have q mn much larger than 1 
the details of A mn are not right and the Gaussian fall-off 
factor in Eq. (JSJ) is not exact, the off-diagonal matrix 
elements with large q mn (large changes of energy in com- 
parison to A) become equivalent to zero anyway. The 
simple Gaussian factor provides the same result. The 
diagonal matrix elements D m with m corresponding to 
E m ^ A are equal to the eigenvalues ui m 3> A of the 
initial H . The RG evolution of the matrix elements with 
qmn ~ 1 was calculated numerically. Eq. © is accurate 
for matrix elements with q mn <C 1. The smaller q mn the 
better the accuracy of Eqs. and ©. The key question 
here is how g\ depends on A. 

The best way to explain the behavior of g\ is to use a 
generic example. On the basis of numerical experiments 
with various choices for the parameter b, coupling con- 
stants g\ and h, ultraviolet cutoff A, and infrared cutoff 
Em , one representative case was selected that clearly dis- 
plays key features of the cycle of the coupling constant 
g\ in effective theories. The set of parameters includes 
b = 4, g A = 0, h = tan tt/5, N = 16, and M = -25. The 
value of b was selected as a result of compromise between 
making it large enough for a qualitative analysis (see next 
section) to be sufficiently precise for numerical confirma- 
tion, and keeping it small enough for the entire range 
of A to remain within bounds set by requiring at least 
four digits of numerical accuracy from 4th order Runge- 
Kutta integration procedure for relevant matrix elements 
of H\. gA = as an initial condition was selected to 
keep the calculation as simple as possible. Small values 
of gA will be discussed later in connection with asymp- 
totic freedom. The selected parameters are similar to 
those used in Refs. [35|, |36| that analyzed the case with 
h = (asymptotic freedom), to facilitate comparison. 

Once the initial Hamiltonian is chosen, the similarity 
RG transformation is fully determined by Eq. J5J. The 
running coupling constant is defined from the full solu- 
tion through 

gx = 1 — H\mm I Em , (7) 



which is analogous to the definition of the fine structure 
constant in the Thomson limit. The continuous line in 
Fig. n displays the resulting g\. The broken line results 
from an approximate formula to be discussed in the next 
section. 
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FIG. 1: The solid line displays the cycle. The broken line is 
an approximate result explained in Section ITTT1 



It is clear from Fig. ^ that the effective coupling 
constant oscillates with period r = b p , where p = 
7t/ arctan [h) = 5. For example, one can read the hor- 
izontal distances between successive points on Fig. Q] 
that correspond to g\ = 0.5 and divide them by 5 In 6. 
From the numbers that were used to draw Fig. ^ one 
obtains (in the order from the largest to lowest value of 
In A in the figure): 1.0326, 1.0032, 0.9996, 0.9999, 0.9995, 
0.9998, 0.9994, 0.9998. 

It is also clear that the pattern of change of effective 
theories in a single cycle has two characteristic rates. 
Namely, the coupling constant changes gradually over 
almost the entire cycle and then suddenly jumps back 
to its initial value. This is a generic feature. It was 
observed numerically for all combinations of the param- 
eters 1.1 < b < 100, 2 < p < 42, and coupling constants 
<?A < 1, with no indication of disappearing outside this 
range. For example, if p — 15, i.e., for h — tan(7r/15), 
the cycles are 3 times longer than in Fig. ^ but the 
jump takes about the same distance along the horizontal 
axis. For large p that considerably exceeds N—M+l and 
corresponds to an extremely small skew symmetric imag- 
inary part in the interaction Hamiltonian, the cycle is so 
long that one can hardly distinguish the cyclic behavior 
of the coupling constant in the model with N = 16 from 
asymptotic freedom (see Fig. [3] in the next section). The 
next section discusses the mechanism of the cyclic change 
and shows how the cases of limit cycle and asymptotic 
freedom are connected. 
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III. THE CYCLIC CHANGE OF g x 

The pattern of incremental changes followed by a sud- 
den big jump that reverses the result of the gradual 
changes and starts a new period of incremental variation, 
is a characteristic feature of the model. The pattern is 
of interest by itself and may provide useful analogies in 
studies of other systems that display such behavior. 

Matrix elements of the Hamiltonian H\ can be written 
in the spectral form, 

JV 

Hxmn = ^2 WfcVfcm(A)l/'fc n (A) , (8) 
k=M 

where cuk is the fc-th eigenvalue of the model and i\)km 
denotes m-th component of the corresponding effective 
wave function. The cigenstates are numbered from M 
to N in the order in which their eigenvalues appear on 
the diagonal when the similarity procedure is carried out 
all the way down to A = 0. The wave functions are 
normalized by the condition J2m=M l^fcmWI 2 = 1. The 
eigenvalues are independent of A. 

The spectral decomposition implies through Eq. Q 
that 

1 N 

SA = ^^(X) , (9) 

M k=M 

where afc(A) = \tp k M(X)\ 2 is the probability of the low- 
est energy component in the eigenvector number k in 
the effective basis corresponding to A. Taking into ac- 
count that in the flow of H\ towards small A large eigen- 
values appear one after another on the diagonal and 
off-diagonal matrix elements go to zero, one can con- 
clude that the numbers a^(A) are reduced to zero one 
after another starting from large k. The eigenvector 
components are expected to vary with A as the inter- 
action does, and the approximate Eq. |J5j suggests that 
afc(A) ~ a fc (oo) exp [—2(E). — E M ) 2 /\ 2 ]. A simple fur- 
ther consideration provides some insight concerning more 
precise description of the behavior of a^(A) as a function 
of A. 

Namely, in the case of N = M + 1, one can solve Eq. 
|J5J analytically and the solution says that the small com- 
ponents of the eigenvectors indeed fall off as a Gaussian 
function, but the rate is given by exp [— (u>k — lom) 2 / X 2 ], 
i.e. the eigenvalues of Ho are replaced by the eigenvalues 
of H . Thus, knowing eigenvectors and eigenvalues of H 
one can write 

a fe (A) = c fe (A) a fc (oo) exp [-2(u k - uj m ) 2 /X 2 } . (10) 

Cfc(A) is an associated normalization factor, ft is com- 
puted using the condition that the similarity transfor- 
mation does not change the norm and assuming that the 
RG evolution of the effective wave functions is given by 
the simple formula 

ipkm(X) = c fc (A) ij; km (oo) exp [-(oJk - w m ) 2 /A 2 ] . (11) 



Phases of the complex wave-function components 
ipkm(oo) increase in steps with m for m > k with a period 
equal to the number of the steps needed to obtain 2ir. A 
similar pattern occurs in the three-body dynamics |l8j . 

Eq. © with formula (|llfl for ipkm{X) is now used to ap- 
proximate the entire RG evolution of the effective Hamil- 
tonian matrices obtained from solutions of Eq. J3J): the 
effective wave functions evolve with A through the ex- 
ponential factors exp [— (uok — uj m ) 2 /X 2 ] and the normal- 
ization coefficients Cfc(A). The fall-off pattern of the ap- 
proximate Eq. l|l(Jf) is compared with exact numerically 
computed behavior of the wave functions in Fig. fusing 
the example of 5 states forming a cycle. 
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FIG. 2: Full lines show the A-dependence of lowest energy 
components, i.e., &\ — &k(\)/ak(oo), for 5 eigenvectors be- 
longing to one cycle with k — —6, —5, —4, —3, and —2, from 
left to right, respectively, and the broken lines display results 
of the approximate Eq. 

A comment is required concerning the order of num- 
bering of eigenstates, which is related to why 6 = 4 was 
chosen in the example. The need for a comment stems 
from the intrinsic complexity of the eigenvalue problem 
and RG equations. A natural ordering of eigenvectors is 
found by tracing the sequence in which eigenvalues ap- 
pear on the diagonal when A is lowered. It turns out that 
the ordering depends on the size of b when the other 
parameters in the model arc kept fixed. For example, 
there is a change in the ordering of eigenvalues between 
b = 2.15 and 2.20 in the case considered here. The 
entire Hamiltonian dictates where on the energy scale 
the negative bound-state eigenvalues are to appear. In 
other words, there is one negative eigenvalue per each 
and every quartet of successive positive eigenvalues but 
the modulus of the negative eigenvalue can wander across 
one of the positive eigenvalues on the diagonal when b is 
changed, b — 4 was chosen to stay clear of this feature 
near b ~ 2.2. But the interference among the evolv- 
ing wave functions is not avoided entirely because the 
bound state wave functions spread coherently over the 
basis states and have more than one significant compo- 
nent. The interference is the origin of the difference be- 
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tween the qualitatively correct Eq. i|10|l and the actual 
dependence of wave functions on A, visible in Fig. [3 In 
any case, all eigenvalues occur in sequences of powers of 
r = b p . There are p — 1 such geometric sequences of 
positive eigenvalues and one such geometric sequence of 
negative eigenvalues. The sequence of negative eigenval- 
ues is shifted on the logarithmic scale with respect to the 
sequences of the positive eigenvalues. Moduli of the neg- 
ative eigenvalues are close to geometric averages of the 
neighboring positive eigenvalues. It is enough to find nu- 
merically where one of the moduli of negative eigenvalues 
fits in among the positive eigenvalues to locate moduli of 
all other negative eigenvalues on the energy scale (a nu- 
merical pattern loses accuracy near the cutoffs M and 
TV). 

Using Eq. © with input from Eq. (|10|l . one obtains 
the broken line shown in Fig. 2] There is only a slight 
discrepancy between the broken line and the actual cy- 
cle, which follows from the complex interference pattern 
that is not fully captured. An additional benefit of the 
approximate formula is that it allows one to study be- 
havior of the coupling constant for large TV, M|, p, and 
6 approaching 1 with much less numerical effort than the 
Wegner equation itself. Although it is easy to under- 
flow or overflow computer accuracy in studies of the log- 
arithmic effects, some preliminary studies indicated that 
the effective coupling constant g\ makes its rapid but 
smooth transition at the scale of binding for large p and 
b approaching 1 without any significant increase in size. 
The maximal value of the coupling constant may even 
decrease when b is lowered toward 1. The width of the 
transition region always appears to match the width of 
the bound-state wave function. 

It is clear from Fig. Q]that the basic features of the cy- 
cle in behavior of g\ as a function of A can be described 
in the following way. The similarity RG flow eliminates 
interactions that change energy by more than A and this 
forces eigenvectors to have their significant components 
squeezed on the H energy scale to the region of size 
A around the corresponding eigenvalues. The eigenval- 
ues appear on the diagonal of H\ one after another and 
the eigenvectors end up contributing only to the diagonal 
matrix elements, being eliminated from further evolution 
entirely. On the other hand, the elimination of p states of 
a full cycle brings the low-energy matrix elements, and 
therefore also the coupling constant, to the same val- 
ues. A similar effect of elimination of states was shown 
already in Refs. [lEl 

using the Gauss elimination pro- 
cedure as the RG transformation. In fact, elimination of 
p high-energy states produces a cyclic coupling constant 
behavior because 

71 

w fc a fc (oo) = 0. (12) 

k—n—p+l 

This identity is satisfied in the initial model Hamilto- 
nian for M <C n <C TV with great accuracy independently 
of the RG procedure. The same pattern matters when 



similarity changes afc(oo) to a/^A) and low energy com- 
ponents of high-energy states are eliminated one after 
another by the Gaussian weight that is characteristic of 
Wegner's equation. Since there is only one bound state 
per cycle, one term in the sum of Eq. 1)12(1 is a negative of 
all the others in a cycle together. This is why the incre- 
mental change in the coupling constant induced by RG 
evolution of p — 1 positive-eigenvalue eigenstates is un- 
done by evolution of one negative-eigenvalue eigenstate. 

The actual width of the rapid transition of g\ as a func- 
tion of A must be related to the width of the bound-state 
wave function on the energy scale (in quantum field the- 
ory, the scale of eigenvalues of Hq corresponds to the scale 
of kinematic momentum variables). For large 6, only a 
few components of the bound-state wave function have 
significant size and the transition occurs over a few pow- 
ers of b. When b is lowered, more components matter. 
The width of the bound-state wave function must deter- 
mine the size of the rapid transition region because all 
significant components of the bound-state must be elim- 
inated in order to produce the bound-state eigenvalue 
on the diagonal. Since a simultaneous rescaling of the 
parameter A and all matrix elements of a cyclic Hamilto- 
nian by r recovers the same array of matrix elements, the 
widths of successive bound states in the cycle also change 
like successive powers of r when A is lowered. Thus, the 
width of the bound states that are eliminated when A 
passes the scale of binding is constant on the logarithmic 
scale. The cyclic variation of g\ must therefore occur 
with a constant width of the rapid-fall-off region on the 
logarithmic scale. 

When p becomes very large, the incremental steps be- 
come very small in size in comparison to the whole cycle 
(for fixed b). But after the coupling g\ grows to values 
much larger than h, its further change is governed by the 
positive spectrum practically independently of h. The 
same behavior of g\ for A near the scale of binding is 
obtained as in the case of asymptotic freedom in Refs. 
HE HH> where h = 0. Fig. |3] illustrates how g x be- 
haves in the two cases, both specified by b = 4, TV = 16, 
M = -25. 

The first case corresponds to a cycle with h = 
tan (n/50) and <?a = 0. In this case, there happens to be 
just one bound state with energy E = —7.644479 10 -6 
(if TV is increased to 66 and all other parameters are 
the same, another bound state with eigenvalue E' = 
—9.690529 10 24 is created in the first case and E' — 
4 50 i?). The second case corresponds to asymptotic free- 
dom. The second case also has just one bound state with 
the same binding energy but h — and <?a = 0.04000228. 
In the second case, the bound-state is present due to a 
non-zero c/a- It is visible in Fig. [3] that as soon as g\ 
evolves toward decreasing values of A in the first case 
(limit cycle, solid curve) to values much larger than h, 
its further RG evolution toward still smaller values of A 
is governed by the size of g\ as in the second case (asymp- 
totic freedom, dashed curve). 

The rapid transition in the dependence of g\ on A when 
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A passes through the scale of binding corresponds to the 
well-known infrared singularity in QCD (infrared slavery) 
that occurs when the momentum scale approaches Aqc d 
in perturbation theory. In the non-perturbative calcula- 
tion in the model, the singularity appears smoothed be- 
cause it is spread over the width of a bound-state wave 
function. 

The asymptotically free model is analogous to a cycle 
of a very large range (very large p for fixed b, which 
means very small h) when g& is set to a finite value in 
the initial Hamiltonian so that the RG evolution brackets 
the sudden-jump area on the energy scale due to a bound 
state. If g\ ^ in the cyclic model, the spectrum is 
slightly changed. It is shifted on the logarithmic scale and 
the range of the cyclic variation of the coupling constant 
is changed. The qualitative features for scales near the 
range of binding remain the same. 

The key difference between the cycle and asymptotic 
freedom shows up when the coupling constant g\ ap- 
proaches 0. In the case of asymptotic freedom, the cou- 
pling tends to as a fixed point when A increases toward 
infinity. In the case of a cycle, when g\ drops in size 
below h, the latter takes over in the dynamics and the 
evolution toward larger A is altered. 
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FIG. 3: Behavior of g\ near the scale of binding in two 
cases: limit cycle and asymptotic freedom. The solid curve is 
obtained for the cycle with p = 50 and the dashed curve for 
asymptotic freedom with g& ~ 0.04 (see text for details). 



IV. DISCUSSION 

Limit cycles of realistic effective quantum theories that 
originate from models with divergences are difficult to 
study because of their complexity. The complexity can 
also obscure an underlying RG structure that may be 
close to a limit cycle, as occurs in a number of cases 
mentioned in the Introduction. The model discussed 
here is so simple that the cycle of effective theories in it 
can be exposed and analyzed without great effort. The 
study shows in a few steps how the similarity RG pro- 
cedure may be set up for investigating limit cycles of ef- 
fective theories. The cyclic behavior of the coupling con- 



stant g\ is found by solving Eq. J5J without calculating 
eigenvalues of the model with large cutoffs. The struc- 
ture of counterterms is found using perturbative similar- 
ity RG transformation order-by-order. The transforma- 
tion guarantees that eigenvalues of effective Hamiltonians 
with finite width A are independent of A. 

The same similarity RG procedure that works in the 
case of a limit cycle also works in the case of asymptotic 
freedom. The procedure shows how the two cases are 
related. It also produces an approximate structure of 
solutions of effective theories in the form of Eq. 
This approximate structure works in both cases and thus 
appears to be more generally valid than the models used 
here to identify it. The meaning of this result can be 
described in the following way. 

Suppose a Hamiltonian for some physical system is pro- 
posed. Suppose that the proposed Hamiltonian renders 
results that are sensitive to arbitrary cutoffs that are not 
an intrinsic part of the proposal but a technical element of 
necessity: the proposed interactions spread states all over 
the space considered in the theory up to the cutoffs and 
the model without cutoffs, strictly speaking, does not ex- 
ist. Matrix elements of the model Hamiltonian can be re- 
interpreted according to Eqs. (J5J and 10 as correspond- 
ing to an effective theory. The re-interpretation consists 
of associating the basis states with energies known from 
experiments and multiplying the model wave functions 
for any given cutoff with a form factor analogous to the 
Gaussian factor present in Eq. The form factor 

contains A. For some finite value of A on the order of 
observables of interest, the form factors in the new (ef- 
fective) Hamiltonian and in the wave functions eliminate 
sensitivity to the arbitrary technical cutoffs. The effec- 
tive theory needs adjustment of constants at some suit- 
able value of A. 

In fact, a similar sequence of steps characterizes con- 
struction of effective theories in atomic, nuclear, or parti- 
cle physics in which one introduces form factors to smear 
singular interactions. (The similarity RG evolution of g\ 
in Hamiltonians of this article can be precisely related 
to the evolution of coupling constants in Hamiltonians of 
quantum field theories with form factors, e.g., see Ref. 
[33.) All that the similarity RG procedure provides here 
is a new context of looking at the construction of effec- 
tive theories, especially in the case of limit cycles. In this 
context, when the energy changes described by the effec- 
tive Hamiltonians increase, new bound states of smaller 
sizes (the size corresponds to 1/A) become a part of the 
effective dynamics. Of course, a bound state of a very 
small size cannot be resolved in processes limited to small 
changes of energy. There, it can only play the role of a 
new particle. 

The close connection between two models: one with 
a limit cycle of a large scale period r = exp(7r//i) and 
another one with asymptotic freedom, exposed here by 
the similarity RG procedure, provides further insights 
concerning a very large range of scales. It suggests that 
Hamiltonians with extremely small imaginary matrix el- 
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ements (extremely small coupling constants analogous to 
h after In b is factored out) can produce a new genera- 
tion of binding effects where asymptotic freedom alone 
predicts no new physics. The limit cycle structure thus 
appears to suggest a scenario in which a dynamical corre- 
lation can emerge across a very large range of scales due 
to an extremely small coupling constant in a well-defined 
theory; the smaller the coupling constant h, the larger 



the scale period r. This is also a good reason for under- 
taking careful studies of limit cycles of effective theories 
using the similarity RG procedure and other means. 
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